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Abstract 



Oh' In this paper we suggest gauge invariant discretization of Poincare 



quantum gravity. We generalize Regge calculus to the case of Riemann- 
Cartan space. The basic element of the constructed discretization is 



^vq I piecewize linear Riemann-Cartan space with flat pieces of hypercubic 

form. We consider the model with squared curvature action and cal- 
culate the correspondent lattice action. We construct local measure 
over the dynamical variables of the lattice model. 
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1 Introduction 

In order to put quantum gravity theory on the lattice one should construct 
such discretization of the continuum model, which preserves as much symme- 
try of the original continuum model as possible. That's why Regge quantum 
gravity is considered to be one of the most natural lattice realizations of quan- 
tum gravity. The main success of Regge calculus is that it is manifestly gauge 
invariant, i.e. it preserves the whole symmetry of the original model. This 
is due to the fact, that in this lattice theory piecewize - linear Riemannian 
manyfold plays the role of the dynamical variable. (In the original contin- 
uum theory Riemannian manyfold of any possible form plays the role of the 
dynamical variable.) In other words, we approximate any given Riemannian 
manifold by piecewize - linear Riemannian manyfolds. This discretization is 
gauge invariant by the construction. 

In order to define measure over link lengths in Regge quantum gravity one 
may start from the invariant metric on the space of continuum Riemannian 
geometries [U El E] : 

\\S9r = lJ d''^\/\9\i9'"'9'"' + 9'" 9"' + Cg^^g''^)5g,,5g,, (1) 

where C is the arbitrary constant such that C ^ —2/D. Then we may 
postulate that the correct measure on the space of Riemannian geometries 
corresponds to this metric in the same sense as the measure dx on R^ cor- 
responds to the metric ||5xp = ((5x)^. Unfortunately, the shift from finite 
dimensional case to infinite dimensional case is not clear. Namely, we may 
generalize the finite dimensional formula for the measure correspondent to 
the given metric. Thus we obtain the continuum measure correspondent to 
the metric ([T]) (for more details see, for example, |5| and references therein): 



Dg = Ii^{^\g{x)\rn,>,dg^,{x) (2) 

Here the choice a = {D — 4){D + 1)/4 corresponds to the super-metric tensor 
in its form chosen in [Hll] while the choice a = —{D + 1) was considered in |2]. 
Obviously this ambiguity is related to the fact, that the product over points 
IIj; is not well defined if x belongs to R^. Actually, the precise meaning could 
be given to ([2]) only when the discretization of Riemannian space is chosen. 
In some of the papers devoted to Regge quantum gravity, the following 
measure over link lengths was considered (see, for example, [5j): 

Di = [u^v:m,dii]e{k,), (3) 



where the product Ux is over the simphces while the product Uij is over the 
vertices of the simphcial manyfold. Q^kj) is the step function, which provides 
the triangle inequalities at each simplex. Vx is the volume of the simplex x, 
and lij is the length of the link that connects vertices i and j. 

According to [5] the derivation of the measure ([3]) is as follows. At each 
simplex V^Ilijdlfj{x) is the direct discretization of the expression 



\g{x)\yil^>,ydgfj,iy{x) (where i and j are the vertices of the given simplex, 
and the constant factor is omitted). In the system of simplices glued to- 
gether the link lengths are not independent. Therefore, the correspondent 
constraint should be imposed. In [5] this constraint was chosen in the form: 
IlxyIlijS{lfj{x) — lfj{y)), where the product II^-j^ is over the pairs of neighbor- 
ing simplices while the product Uij is over the pairs of vertices of the side, 
which is common for x and y. Thus we obtain the lattice measure in the 
form ([3]). 

We mention here again, that the derivation of the continuum formula ([2]) 
and even its own form cannot be considered as completely mathematically 
rigorous. There is also an ambiguity in the choice of the constraint, which 
is imposed on the link lengths in the system of simplices glued together. 
Namely, such a constraint may be choosen in the form Q[lij]llxyllij6{lfj{x) — 
If Ay)), where ^\lij\ depends upon link lengths. A priory the particular form 
of such functional is not known. Therefore, although the choice (jS]) seems 
to be quite natural, it is not derived rigorously from the expression ([1]) for 
the metric on space of continuum Riemannian geometries. This results, in 
particular, in the ambiguity in the choice of a. 

A possible solution for this problem would appear if there exists a symme- 
try in the continuum theory, which (after transferring to discretized model) 
fixes the form of lattice measure. Unfortunately, the existence of such a sym- 
metry within Regge quantum gravity is not clear. Therefore, in the present 
paper we suggest to consider its generalization based on the so-called Poincare 
quantum gravity. We hope that the mentioned problem can be solved within 
this generalized model. This is the theory, in which the dynamical variable is 
the Riemann - Cartan manyfold. (Riemannian manyfold is the limiting case 
of Riemann - Cartan manyfold with vanishing torsion.) 

It is worth mentioning, that the topics related to the discretization of 
Poincare gravity were already considered. Say, in P [7] the independent 
5*0(4) connection was introduced into the Regge calculus. The connection 
is singular and lives on the sides of the simplices. In principle this con- 



struction resembles the one suggested in the present paper. However, here 
we use hypercubic lattice, which is much more useful, than the simplicial 
one. Moreover, in the present paper we construct simple natural Poincare 
invariant measure over the dynamical variables, which was not done in ^ [7] . 
We also mention, that in [9], [7] there was no consideration of the squared 
curvature terms in the action. 

Poincare gravity on the hypercubic lattice was considered, say, in [HI [H 
[Tot [TT| [T2| [13] . In these papers lattice discretization was constructed in 
the conventional way via direct discretization of the theory. Therefore, the 
correspondent constructions were not gauge invariant (in the gravity models 
the gauge group is the group of general coordinate transformations)!!!. In 
the mentioned papers some of the considered lattice theories were invariant 
under lattice gauge transformations, which belong to the Poincare group. 
However, the existence of this invariance does not mean that the constructed 
discretizations are invariant under general coordinate transformations. For a 
brief comparison of our approach with the approaches of [8l HI [101 El Il2l [13] 
see section 3 of the present paper. 

In the present paper we generalize Regge construction to Riemann-Cartan 
spaces. In our approach we approximate any given manyfold by a piecewize- 
linear Riemann-Cartan manyfold (with hypercubic elements). The crucial 
difference between our approach and the approachs of [H [9], [TOj [HI [121 [13] is 
that we use geometrical construction. I.e. we discretize the original many- 
fold by the invariant objects (with respect to the group of general coordinate 
transformations). As a result our lattice model is gauge invariant by con- 
struction. 



2 The discretization of Riemann-Cartan space 

Thus, we consider Riemann-Cartan space, in which both 5*0(4) connectioiu 
and the inverse vierbein are the dynamical variables. We do not require van- 
ishing of torsion (which would lead to appearance of Riemannian manyfold) 
or vanishing of curvature (which would lead to appearance of Weitzenbock 



^It was shown, for example, that in ordinary gauge theories the discretization, which is 
not gauge invariant, is not appropriate (say, in the correspondent lattice nonabclian gauge 
models there is no confinement of fundamental charges [l4]). 

^Further we imply, that Wick rotation to Euclidean signature is performed. So, we 
deal with Euclidean path integral formalism. 



space). The discretized space is in itself a Rieinann-Cartan space. It is com- 
posed of flat pieces connected together. We consider the case, when each such 
piece has the hypercubic form. Further we refer to hypercubes as to elements 
of the lattice. Form of the lattice elements is fixed by the set of vectors e^ 
that connect the center of the element with its vertices. The expression of 
e^ through elements of the orthonormal frame f^ {A = 1,2,3,4) (common 
for all lattice elements) is one of the basic variables of the construction. So 
we have 



e 



.-E<fA (4) 

A 

(Everywhere space - time indices are denoted by Greek letters contrary to the 
tetrad ones.) We imply here that vectors e^(/i = 0, 1, 2, 3, 4) are independent. 
The other vectors e^(/i = 5, ..., 15) are defined in such a way, that opposite 
sides of the lattice element are parallel to each other. The hypercubic lattice 
is periodic and the position of the starting point of each lattice element is 
always denoted by Bq. Vectors e^(/i = 1, 2, 3, 4) point to its neighbors. 

Variables E^ represent translations from the center of the lattice element 
to its vertices. Metric (or vierbein) is implied to be constant inside each lat- 
tice element. Here shift of the center of lattice element by a vector v"^ causes 
transformation of basic variables: E^ -^ E^ + v^, which could be treated as 
gauge transformation with respect to the translational gauge group. It rep- 
resents the translation of the given lattice element within the correspondent 
local map. 

In addition to the translational connection, which is defined by the set of 
variables E^, each shift from one lattice element to another is accompanied 
by the rotation in the four - dimensional tangent space. In other words, there 
is the SO {4) connection, which is singular and is concentrated on the sides 
of lattice elementqj- We denote by Uu the SO (A) matrix, which is attached 
to the side that is common for the lattice elements I and J. 

The constructed Riemann - Cartan space has singular connection. In 
this case definitions of curvature and torsion become ambiguous. Therefore 
we must fix one of the definitions. For the details we refer to the Appendix, 
where one of the definitions is used in order to calculate torsion and curvature 



^In this paper we do not consider parity transformations, correspondent to the matri- 
ces diag{±l, ±1, ±1, ±1} with negative determinant since these transformations do not 
appear in the infinitesimal form of Poincare transformation. However, the considered 
model could be easily generalized in order to include parity transformations if we use 0(4) 
matrices instead of 50(4) ones. 



on the piecewize-linear manyfold. 

The gauge transformations with respect to the whole Poincare group are 
represented by translations and 5*0(4) rotations of the lattice elements, that 
result in the following change of basic variables: 

< - 0B< + v\ (5) 

where is the rotation matrix and v is the vector that represents translation. 

3 The action 

In this paper we consider the following action: 

S = f{a{RABRAB--R^)+f3R^-irnlR + Xm^p}\S\d^x 

+ 5mljTicT^cW\ (6) 

where \£\ = detSf^, Sf^ is the inverse vierbein, the tetrad components of Ricci 
tensor are denoted by Rab, and R is the scalar curvature. Coupling constants 
a, P, 7 and A are dimensionless while rrip is a dimensional parameter. The 
second term is added in order to suppress torsion at 6 -^ oo. So, at (5 ^ oo 
we arrive at the model considered in [15]. This model is renormalizable and 
asymptotic free at certain values of the coupling constants. Thus the model 
with the action iQ is worth considering. More arguments in favor of this 
point of view could be found in [16j. 

The action (jH]) can be calculated on the piecewize - linear Riemann-Cartan 
manyfold composed of flat hypercubic cells. Below we represent the result- 
ing expression. For the details of the definition and calculation of singular 
curvature and torsion on this manyfold see Appendix. 

In order to rewrite in a useful form the expressions (TT4|) and (TTSl) for 
lattice curvature and toirsion we drop to the dual lattice. Then our rotation 
matrices are attached to links while the inverse vierbein is attached to sites. 
Let us denote by Un{x) the matrix correspondent to the link, which begins 
at the site x and points to the direction /x (// = ±4, ±3, ±2, ±1). We denote 
by Q^i,{x) = U^{x)... the product of link matrices along the boundary of the 
plaquette, which is placed in the (/ii^) plane. The inverse vierbein, which is 
attached to the site x, is £^ = E^ — E^, /i = 1, 2, 3, 4. For negative values of 



/i we define £^ = —£a- The inverse matrix for positive values of indices is 
denoted by ^^(x) = {£{x)~^}^. We also expand this definition to negative 
values of indices: ^^(x) = sign(/i){[£^(a;)]~^}j4 . We shall denote by Ax^ the 
shift on the lattice by one step in the /i -th direction (Ax_^ = — Ax^). So, 
X + Ax^ is the site which is obtained via the shift from the site x by one 
lattice spacing in the direction /i while x — Ax^ is obtained by the shift in 
the opposite direction. Thus, U_^{x) = U^^{x — Ax^). Next, we define 
A^£y{x) = U^{x)£y{x + Ax^) — £y{x). Everywhere we imply summation 
over the repeated indices. The summation over space - time indices /i, i/, ... 
is implied over ±4, ±3, ±2, ±1. 

We introduce lattice tetrad components of torsion and curvature attached 
to the sites of the lattice and to positive or negative directions (in this para- 
graph and in the next paragraph the summation over Greek indices is not 
implied) : 

, ^A 



[RnFBi^) = l^nAFBi^) 

[R'^l(x) = [R^l^^(x) 
R(x) = 5: [R'^l(x) 

11,1/ 

V 



(7) 



Now we are ready to express the action ^ on the piecewise-linear many- 
fold through its parameters E^{x) and U^{x). In order to do this we use the 
expressions for torsion and curvature flT^ and flTSj) . We omit the intermedi- 
ate steps and represent the final result using the notations for lattice torsion 
and curvature introduced above. 

s = i: \£{x)\{j:MRnFBi^mnFBi^) + i^-l<^Mn'ix) 



xGsites Mi^ 



-7mJ[Rn(x)} + SmlY.[T^^]Bci^)[T^TBci^) + ^^p)> (§) 

where \£{x)\ = |det£^(x)| is the volume of the lattice element correspondent 
to the site x. 

Here we introduce lattice couplings a, [3,6, A, 7 that differ from the orig- 
inal ones by the factors, which are formally infinite and come from delta - 



functions in expressions (IT^ and (IT^ for torsion and curvature. We assume 
here, that a certain regularization is made, which makes these factors finite. 
Our supposition is that after the renormahzation each physical quantity may 
be expressed through physical couplings, which differ from the bare ones 
(both lattice and continuum), and the infinity encountered here is absorbed 
into the renormahzation factors. 

In our lattice model the translational connection E^{x) attached to sites 
and the SO (4) connection U^j_{x) attached to links are the dynamical vari- 
ables. The action of the model is expressed in a compact way through these 
variables. It is easy to understand, that (jSj) is the symmetry of the action. 
So, we have lattice model with the direct manifestation of Poincare gauge 
invariance. 

As we have already noticed in the introduction, lattice realization of 
Poincare quantum gravity has already been considered in several papers (see, 
for example, [H El HOI HU Il2l [13] and references therein) . Now we are ready to 
compare our approach with the approaches of the mentioned papers. First, 
in all of them the continuum models were transferred to lattice via appli- 
cation of naive discretization procedure. This means, that certain variables 
attached to the sites and/or links of the hypercubic lattice were considered, 
and it was set up the correspondence between them and the Poincare group 
connection of the continuum theory. In [9l [121 [13] both vierbein and 5*0(4) 
connection on the lattice were attached to links. In [TU] Poincare group was 
considered as a limiting case of de Sitter group, and link variable belongs 
to 5*0(5). In [H [11] the 50(4) connection was attached to links while the 
vierbein was attached to sites of the lattice. The definition of the lattice 
model in [11] is in principle close to the resulting definition of our model on 
the dual lattice. However, the resulting models are not identical. First of 
all, in [TT] the action was of the Einstein - Hilbert form in Palatini formu- 
lation, i.e. it contains only the first power in lattice curvature. The term 
—^rn'^Y.x&sites'^{.^)\^{.^)\ of (IHD would coincide with the action of [11] if in 
the latter the symmetrization over different directions at each site is per- 
formed. In [S] the squared curvature action of general type was considered. 
Both the link 0(4) connection and the site inverse vierbein were considered 
in 4 X 4 spinor representation. The final expression for the lattice action is 
rather complicated and does not coincide with ([8]). 

The crucial difference between the mentioned approaches and the ap- 
proach of the present paper is that we use the regular procedure and ap- 
proximate Riemann-Cartan manyfold via piecewize - linear Riemann-Cartan 

8 



manyfolds. This procedure gives us the discretization of the original contin- 
uum model, which is manifestly invariant under general coordinate transfor- 
mations. Contrary to this, the constructions considered in [8| [9| [TO l ITT| [T2 | [T3] 
violate general coordinate invariance and do not give gauge invariant dis- 
cretization of Poincare gravity. 

4 Measure over discretized geometries 

In order to define measure over dynamical variables in our lattice model we 
use an analogy with QCD, the lattice version of which is considered to work 
perfectly. In QCD there are two kinds of fields: 

1. Quarks and leptons. (The correspondent measure over Grassmann 
variables is well defined and unique.) 

2. The gauge field. The correspondent measure on the lattice is unique 
as it is completely defined via symmetry properties: that is the local measure 
invariant under gauge transformations. 

In our case there are two fields: 5*0 (4) connection and the translational 
connection. So, it is natural to use measures, which are invariant under lat- 
tice realization of the gauge transformation. Our choice of measure is the 
measure, which is simultaneously invariant under lattice gauge transforma- 
tions and is local. 

Each piecewise linear manyfold described above is itself a Riemann - Car- 
tan space. Let the given discretization (with varying E and U) be denoted 
as Ai. Then, we consider the set of correspondent independent variables 
{E^(I); Uij}. Gauge transformation corresponds to the shift of each lattice 
element by the vector v^{I) and its rotation Gi G 5*0 (4). This transforma- 
tion acts as {E^iI);U,j} ^ {QiE^il) + viI)y,QiU,jQ^}. 

The locality of lattice measure means the following. The whole measure 
should be represented as 

Dm{E; u) = nin^DE;f (i)ni,jDf/ij, (9) 

Here the product is over the sides of lattice elements and over the links that 
connect centers of lattice elements with their vertices. The measure over link 
matrices Uu is denoted by DUij. The measure over vectors E^{T) is denoted 
by DE^iJ). We call the lattice measure local if inside each lattice element 
DE^^ for the given /x depends upon E^^ only, and DUu for the given I, J 
depends upon Uu only. It is obvious that this requirement together with 



gauge invariance fixes the only choice of DE^ and DUj,: DE^ = Yl^^^dE^, 
while DU is the invariant measure on SO (4). 

We must mention, that another locality principle can be formulated. Say, 
we may thought that the measure is local if DE^ may depend upon E^ 
with v ^ i^ but it may not depend upon E-^ from another lattice element. 
Then gauge invariance does not fix measure precisely. However, we choose 
the more strong requirement that was described above since it gives us an 
opportunity to fix the only local gauge invariant measure. 

5 Metropolis algorithm. 

It is worth mentioning, that in real numerical calculations it would be useful 
to express each SO {A) link matrix (on the dual lattice) as the function of the 
SL(2,C) matrix. The correspondence is given by the conventional spinor 
representation of 5*0 (4) rotations. Then the invariant measure on SL{2,C) 
generates the invariant measure on SO (4). It is much more easy to simulate 
the SL{2,C) field than the ^0(4) field itself. 

Metropolis algorithm for the simulation of our model can be described 
as follows. At each step of the algorithm the given particular link and one 
of its ends are considered. It is formed the proposition of the link SL{2, C) 
matrix and the 4x4 matrix attached to the given endpoint of the linly. 
Then the terms of the action ([H]), which contain torsion and correspond to 
the neighbors of the given point are calculated. Next, we calculate the terms 
of the action, which contain curvature and correspond to the points of the 
"butterfly" correspondent to the given link. The "butterfly" is the figure 
that consists of all plaquettes with the given link as one of their sides. Let 
the sum of these terms be denoted as S^™. Then let the same sum that was 
calculated using the old values of the proposed variables, be denoted as 5'°^'^. 
The proposition is accepted if 5'°'^™ < 5'°^'^. Otherwise it is accepted with the 
probability exp (5'°''^ — S'"'^™). Thereafter we choose another link and one of 
its ends, and the procedure is repeated. 

In order to accelerate the numerical simulation one can save in the com- 
puter memory the values of all plaquette variables fl^^, the values of all 
lattice derivatives A^S^, and the values of the inverse matrices {^^}^^. This 



^We may fix the gauge and choose Eq — 0. The correspondent Faddeev-Popov deter- 
minant is constant. 
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would lead to an essential economy of CPU time but requires an additional 
memory size. 

6 Conclusions 

In this paper we construct the hypercubic discretization of Poincare quan- 
tum gravity. The main achievement of our construction, which differs it 
from the other discretizations of Poincare gravity, is that this construction 
is manyfestly gauge invariant (with respect to the group of general coordi- 
nate transformations). Namely, we approximate continuum Riemann-Cartan 
manyfold by piecewise-linear Riemann-Cartan manyfolds. The resulting lat- 
tice model deals with the invariant geometrical properties of these piecewise- 
linear manyfolds. This construction is analogous to Regge discretization of 
Riemannian manyfolds and is its direct generalization. 

We consider piecewise - linear manyfolds composed of fiat pieces of hy- 
percubic form. So, the geometry is defined by the forms and sizes of these 
fiat pieces together with the 5*0(4) rotation matrices attached to their sides. 
Actually these variables are expressed through the translational connection 
(which is constant inside each fiat piece) and the 5*0(4) connection (which 
is singular and is attached to the sides of the hypercubes). Therefore, these 
translational and rotational connections are the dynamical variables of the 
discretized model. 

The piecewise-linear manyfold has singular torsion and curvature, which 
are calculated directly and expressed through the sizes of the lattice elements 
and the rotation matrices. 

In addition, we calculated the squared curvature action of rather gen- 
eral form on the piecewise-linear manyfold. It is also expressed through the 
dynamical variables of the discretization. 

Next, we point out that the lattice model may be considered as a lattice 
Poincare gauge theory. Namely, shift and rotation of each piece of the given 
piecewise-linear manyfold with respect to the local map cause the transfor- 
mation of the chosen dynamical variables. The lattice action is invariant 
under this transformation. 

We construct local measure over lattice dynamical variables that is in- 
variant under the lattice Poincare transformations. 

Finally, we obtain gauge invariant lattice realization of Poincare quantum 
gravity, in which the translational connection (4x4 matrix) attached to the 
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sites of the dual lattice and the 5'0(4) matrix (or, equivalently, SL{2,C) 
matrix) attached to links of the dual lattice play the role of the dynamical 
variables. The action of the lattice model contains only the first lattice 
derivatives of the dynamical variables and has a compact and rather simple 
form. Thus the constructed lattice model is expected to be useful for the 
numerical simulations. The correspondent Metropolis algorithm is briefly 
described. 

This work was partly supported by RFBR grants 06-02-16309, 05-02- 
16306, and 04-02-16079, by Federal Program of the Russian Ministry of In- 
dustry, Science and Technology No 40.052.1.1.1112. 

7 Appendix 

Here we represent the expressions for curvature and torsion on the piecewise- 
linear manyfold. 

Connection is singular on the sides of lattice elements. 5*0(4) curvature 
is concentrated on the boneqj. We choose the following integral equation as 
a definition of 5*0(4) curvature. 

2 / n{z, y)R^,{y)n+{z, y)dy^ A dy" = Pexp( f u^dx^) 
-[Pexp( r tUf.dx'')]'^ at |S| ^ (10) 



Here cj^ is 50(4) connectioiij. S is a small surface, that crosses the 
given bone, and |S| is its area. dT, is the boundary of S. Its orientation 
corresponds to orientation of S. Q{z,y) = Pexp{J^ uj^^dx^^) is the parallel 
transporter along the path that connects a fixed point on (9S with the point 
y. We choose this path in such a way, that it is winding around the given 
bone in the same direction as 9S and has minimal length. The integral in 
the right hand side is over the path 9S, which begins and ends at the point 
z. 

It is worth mentioning, that the given definition does not contradict with 
the conventional one in case of smooth connection. And it gives us the possi- 



^On the hypercubic lattice we refer to the plaqiiettes as to the bones m order to make 
an analogy with the siniplicial Regge calculus. 

^Uij — Pexp(J w^da:'^), where integral is over the path of minimal length, which 
connects centers of lattice elements / and J. 
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bility to calculate curvature in the case of the constructed singular piecewise 
- linear manyfold. 

Let us fix the given lattice element. Inside it lattice curvature is equal to 

(11) 

Here the sum is over the bones that belong to the given lattice element. The 
integral is over the surface of the bone. The product UiUj j. of the rotation 
matrices is along the closed path around the given bone b, which consists 
of links that connect centers of the lattice elements. Here we imply, that 
this closed path begins within the given lattice element and has the minimal 
lattice length. 

Now let us calculate torsion, which is concentrated on the sides of lattice 
elements. The torsion field T^ is defined by the integral equation 

ni{z,y)T;^^{y)dy''Ady''= f ni{z,y)bl{y)dy^ (12) 

Here hf^{x) is the field of inverse vierbein, which is expressed through our 
variables E^ inside each lattice element if the given parametrization of the 
lattice element is chosen. 

This equation is satisfied with the following expression (which is valid 
within the lattice element I): 

K(V) = El/" """-'""^i^'^-'-' dx- A dx' A dx'S^^\y - X) 

sesides •^^'^^ ^■ 

Jx&s D\ 

Here the first integral in the sum is over the given side s seing from the 
neighbor lattice element J^ (the side s is common for I and J^). (We imply 
that in (1131) the given lattice element and all its neighbors have the common 
parametrization. ) 

Let us define inside each lattice element the following variables: £f^ = 
E^ — Eq, /i = 1, 2, 3, 4. Also we denote by S'^ elements of the inverse matrix 
S^^ In tetrad components we have: 

R^FBiy) = ^ii^F E / e,.,.dx^Adx'^6^'\y-x) ^"^^''^'^^ ~ ^"^^''''"^^^^'^ . 

(14) 
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Torsion is expressed as 



sGstdes 



Here £i^ is calculated inside the given lattice element I while S^^^ is calculated 
within its neighbor J^. 
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